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Abstract 

Collaborative beamforming (CB) is a power efficient method for data communications in wireless 
sensor networks (WSNs) which aims at increasing the transmission range in the network by radiating 
the power from a cluster of sensor nodes in the directions of the intended base station(s) or access 
point(s) (BSs/APs). The CB average beampattern expresses a deterministic behavior and can be used 
for characterizing/controling the transmission at intended direction(s), since the mainlobe of the CB 
beampattern is independent on the particular random node locations. However, the CB for a cluster 
formed by a limited number of collaborative nodes results in a sample beampattern with sidelobes 
that severely depend on the particular node locations. High level sidelobes can cause unacceptable 
interference when they occur at directions of unintended BSs/APs. Therefore, sidelobe control in CB 
has a potential to increase the network capacity and wireless channel availability by decreasing the 
interference. Traditional sidelobe control techniques are proposed for centralized antenna arrays and, 
therefore, are not suitable for WSNs. In this paper, we show that distributed, scalable, and low-complexity 
sidelobe control techniques suitable for CB in WSNs can be developed based on node selection technique 
which make use of the randomness of the node locations. A node selection algorithm with low-rate 
feedback is developed to search over different node combinations. The performance of the proposed 
algorithm is analyzed in terms of the average number of trials required to select the collaborative nodes 
and the resulting interference. Our simulation results approve the theoretical analysis and show that the 
interference is significantly reduced when node selection is used with CB. 
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I. INTRODUCTION 

Wireless sensor networks (WSNs) have become practical technology due to the production of 
low cost, low-power, and small size sensors. Different applications such as habitat and climate 
monitoring, detection of human/vehicular intrusion, and etc. are increasingly employing WSNs 
[[ll. Such applications require sensor nodes to be deployed over a remote area to collect data 
from the surrounding environment and communicate it to far base stations or access points 
(BSs/APs). As a result, the challenges faced in the WSN applications are quite different from 
that of considered in the applications of the traditional wireless ad-hoc networks [|2l. These 
differences can be summarized as follows. 

(i) Typical WSN is densely deployed and may consist of thousands of sensor nodes. 

(ii) The network geometry changes all the time due to failure of sensor nodes or deployment 
of new sensor nodes. 

(iii) Sensor nodes are battery-powered and the battery often cannot be replaced. Thus, the sensor 
node life time is limited by the battery lifetime. 

(iv) Sensor nodes have simple hardware with limited computational capabilities and small mem- 
ory in order to keep the production cost of the sensor node reasonable. 

(v) Sensor nodes can fail easily. Thus, it is desired that the WSN performance does not depend 
on individual sensor nodes. 

(vi) Data and traffic models in WSNs depend on the application, and usually the data is 
redundant, while the traffic has low-rate burst nature. 

(vii) Sensor nodes in WSNs are usually deployed at the ground level and have no mobility. Thus, 
the channel path loss for individual node is high and the channel variations are slow. 

Practical communication schemes for WSNs should overcome the problem of limited trans- 
mission range of individual sensor nodes, while being distributed and scalable. Moreover, for 
designing such communication schemes, power consumption and implementation complexity 
issues have to be taken into account as the most significant design constraints for WSNs. 

To address the aforementioned issues, the inherent high density deployment of sensor nodes 
has been used to introduce collaborative beamforming (CB) for the uplink communication to a 
BS/AP Q, [lU. CB extends the transmission range of individual sensor nodes by using a cluster 
of sensor nodes in a power-efficient way. Particularly, sensor nodes from a cluster of nodes 
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act collaboratively as distributed antenna array to form a beam toward the direction(s) of the 
intended BS(s)/AP(s). Given that each sensor node is equipped with a single omnidirectional 
antenna and operates in half-duplex mode, CB is performed in two stages. In the first stage, the 
data from source node(s) in a cluster is shared with all other collaborative nodes, while in the 
second stage, this data is transmitted by all sensor nodes simultaneously and coherently. In the 
latter stage, sensor nodes adjust the initial phase of their carriers so that the individual signals 
from different sensor nodes arrive in phase and constructively add at the intended BS/AP. In 
this way, CB is able to increase the area coverage of WSNs and, therefore, can be also viewed 
as an alternative scheme to the multi-hop relay communications. However, as compared to the 
multi-hop relay communications, CB brings the following advantages. 

(i) For CB, there is no dependency of communication quality on individual nodes. Thus, the 
communication link is more reliable. 

(ii) CB distributes the power consumptions over large number of sensor nodes and balances the 
lifetimes of individual nodes O. 

(iii) CB enables to create a direct single-hop uplink to the intended BS(s)/AP(s). Thus, it reduces 
the communication delay and data overhead. 

(iv) CB achieves higher connectivity than that of omnidirectional transmission with the same 
transmit power jU. 

In order to implement CB, the following issues related to the the distributed nature of the 
WSNs have been addressed. Distributed schemes for estimating the initial phases of the local 
node oscillators in WSNs have been introduced in H-lIll. These schemes allow to achieve phase 
synchronization among all collaborative nodes in a cluster of WSN. Moreover, to minimize the 
time required for multiple sources to share the data among all sensor nodes in a cluster, a medium 
access control-physical (MAC-PHY) CB scheme which is based on the medium random access 
has been proposed in ifTOll . 

Although, the above mentioned implementation issues for CB have been positively addressed, 
one more concern is that the random sensor node locations result in a random beampattem 
which depends on the actual locations. The effect of the spatial sensor node distribution on the 
directivity of the CB beampattern has been studied in Q (see also ifTTI ) for the case of uniform 
sensor node distribution and in [|4l|.[[T2|. lfT3ll for the case of Gaussian sensor node distribution. 
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Although it has been shown for both aforementioned node distributions that the CB sample 
beampattem has a deterministic mainlobe which is independent on the random sensor node 
locations, the sidelobes of the CB sample beampattern are totaly random and can be described 
only in statistical terms ||3l, BUl, [[T3l . In addition, the aforementioned multiple access scheme of 
ifTOl . which minimizes the time required for sharing the multiple source data, results in higher 
sidelobes even for the average CB beampattem. All these can lead to high interference levels 
at the directions of unintended BSs/APs. Therefore, the sidelobe control problem arises in the 
context of WSNs. Indeed, achieving a sample CB beampattern with lower sidelobe interference 
at unintended BSs/APs has the potential to increase the WSN throughput [fT4|. 

Due to the inherent distributed nature of WSNs, the sidelobe control should be achieved with 
minimum data overhead and knowledge of the channel information. Unfortunately, traditional 
sidelobe control techniques ifTSl . ifTTl cannot be applied in the context of WSNs due to their 
unacceptably high complexity and the requirement of centralized processing. Indeed, to apply the 
centralized beamforming weight design in the WSNs, a node or BS/AP has to collect the location 
and channel information from all sensor nodes. It can significantly increase the overhead in the 
network and nullify the above mentioned advantages of CB. Note that for the same reasons the 
recently developed network beamforming techniques |fT9l , [|2T1l are restricted to the applications 
in the relay networks only and it is impractical to apply them for WSNs. 

In this papeiQ, we develop a technique for sidelobe control in CB for WSNs which is based 
on sensor node selection. Such technique makes use of the randomness of node locations and is 
distributed and scalable as well as it has low data overhead. Moreover, as compared to the optimal 
beamforming weights assignment, our sidelobe control technique which is based on the phase 
synchronization and simple node selection is more robust to the channel/phase errors. In addition, 
it helps to balance the life times of all sensor nodes since the corresponding beamforming weights 
have the same magnitude for all nodes. For the sidelobe control technique, a node selection 
algorithm with low-rate feedback is developed to search over different node combinations. The 
performance of the proposed algorithm is analyzed in terms of the average number of trials 
required to select the collaborative nodes and the distribution of the resulting interference. 

'Some preliminary results have been also reported in 1181 . 
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The paper is organized as follows. System and signal models are introduced in Section II. 
A new sidelobe control technique for CB in WSNs is developed in Section III where the node 
selection algorithm is also summarized. The performance characteristics of the proposed node 
selection algorithm such as the average number of trials required to select the collaborative 
nodes and the resulting interference are studied in Section IV. Section V reports our simulation 
results and is followed by conclusions in Section VI. Proofs of some results in the paper are 
summarized in Appendixes. 

II. SYSTEM AND SIGNAL MODEL 

A. System Model 

We consider a WSN with nodes randomly placed over a plane as shown in Fig. [T] Multiple 
BSs/APs, denoted as V = {do, di, d2, ■ ■ ■ , do}, are located outside and far apart from the 
coverage area of each individual node at directions (po, Lpi, (^2, ■ ■ ■ , '^d, respectively. Uplink 
transmission is a burst traffic for which the nodes are idle most of the time and have sudden 
transmissions. Thus, we adopt a time-slotted scheme where nodes are allowed to transmit at the 
beginning of each time slot. The downlink transmissions are mostly the control data broadcasted 
over a separate error-free control channels. The BSs/APs can use high power transmission and, 
therefore, the downlink is less challenging and can be organized as direct transmission. 

We assume that due to the limited power of individual nodes, direct transmission to the 
BSs/APs is not feasible and sensor nodes have to employ CB for the uplink transmission. 
The distance between nodes in one cluster of WSN is small so that the power consumed for 
communication among nodes in the cluster can be neglected. Each sensor node is equipped 
with a single antenna used for both transmission and reception. To identify different nodes and 
BSs/APs, each node or BS/AP to has a unique identification (ID) sequence that is included in 
each transmission. 

At each time slot, a set 5 = {sq, Si, S2, . . . , S5} of source nodes is active. Moreover, only K 
source-destination pairs are allowed to communicate. Here K = min{card(iS), card(r')} where 
card(-) denotes the cardinality of a set, and the A;th source-destination pair is denoted as Sk-dk- 
For source node s^, the coverage area is, ideally, a circle with a radius which depends on the 
power allocated for the node-to-node communication. 
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Let Ai'' be a set of nodes in the coverage area of the node s^- Let us, therefore, select the 
source node as a local origin of the coordinate system used to mark the spatial locations of the 
nodes in the coverage area of s^, i.e. Ai''. The rth collaborative node, denoted as c,., r G Ai'', 
has a polar coordinates (pr,ipr)- The Euclidean distance between the collaborative node and 
a point (A, (p) in the same plane is defined as 

= a/A2 + p2 _ 2pj.A COS(0 - tpk) ^ A- Pr COs(0 - Ipr) (1) 

where A^ in the far-field region. 

The array factor for the set of sensor nodes A^*^ in a plane can be defined as 

AFk{(j)) = J2 V^e^^'e-^'^'-^'^) (2) 

where is the transmission power assigned to each node, 6*^ is the initial phase of the rth sensor 
carrier frequency, 6*^.(0) = ^dr{(f)) is the phase delay due to propagation at the point {A,(p), 
and A is the wavelength. Then the far-field beampattern corresponding to the set of sensor nodes 
M'' can be found as 

BP,{<p)^\AFkm' = 

where | ■ p denotes the magnitude of a complex number. 

It is assumed that the symbol duration is very short as compared to the channel coherent 
time. Therefore, the channel variations can be considered to be quasi-static during one symbol 
transmission. Since sensor nodes are located at the ground level, the large-scale fading is the 
dominant factor for the channels between collaborative nodes in Ai'' and BSs/APs. Then the 
channel coefficient for rth collaborative node which serves A;th source-destination pair can be 
modeled as 

hrk = Clrkbrk (4) 

where brk is the attenuation/path loss factor in the channel coefficient due to propagation distance 
and ttrk is a lognormal distributed random variable which represents the flactuation/shadowing 
effect in the channel coefficient, i.e., ark ~ exp{A/'(m, a^)}. Here m and denote respectively 
the mean and variance, of the corresponding Gaussian distribution. Then the mean and variance 



(3) 



January 19, 2010 



DRAFT 



7 



of the lognormally distributed ark can be found as 

"a,, = B{o,t} = e"'+* (5) 
= £{»4} = (e"'-l)(e--^) (6) 

where E{-} stands for the statistical expectation. The attenuation/path loss depends on the 
distance between Cr and dk and the path loss exponent. Assuming that all nodes in A^^ are 
close to each other, the pass losses from the nodes in fA^ to the BS/AP are equal to each other, 
i.e., hrk = bk, r E Ai^ [|22l . Moreover, since all BSs/APs are located far apart from the cluster 
of collaborative nodes, the network can be viewed as homogeneous and the attenuation effects 
of different paths can be assumed approximately equal to each other, i.e., hk = h [l23l . Note 
that even if the attenuation effects for different BSs/APs are different, they can be compensated 
by adjusting the gains of the corresponding receivers or the power/number of the corresponding 
collaborative nodes participating in CB. 

B. CB and Corresponding Signal Model 

Consider a two-step transmission which consists of the information sharing and the actual CB 
steps. Information sharing aims at broadcasting the data from one source node to all other nodes 
in its coverage area. Specifically, in this step, the source node Sk broadcasts the data symbol Zk 
to all nodes in its coverage area A^'^, where the data symbol -z^ e C belongs to a codebook of 
zero mean, unit power, and independent symbols, i.e., E {zk} = 0, \z1\ = 1, and E {zkZn} = 
for n ^ k. 

In the case of multiple source nodes sharing their own data with other nodes in their corre- 
sponding collaborative sets of nodes, a collision can occur. To avoid the collision, orthogonal 
channels in frequency, time, or code can be used. However, such collision avoidance causes 
resource loss and can lead to the network throughput reduction, especially if the number of 
source nodes sharing the data is large. Therefore, collision resolution schemes can be used 
alternatively (see for example [[24|. [flSl ). In this case, the information sharing takes only one 
time slot in a random access fashion. Finally, we assume that the power used for broadcasting 
the data by the source node is high enough so that each collaborative node c^. can successfully 
decode the received symbol from the source node Sk- 
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During the CB step, each collaborative node Cr,r E Ai'', is first synchronized with the initial 
phase 9^. = —^pr cos{(fk — ^r) using the knowledge of the node locations (see the closed-loop 
scenario in O). Alternatively the synchronization can be performed without any knowledge 
of the node locations (see [|7]|-[|9l|). For example, the synchronization algorithm of [|7]| uses a 
simple 1-bit feedback iterations, while the methods of ^ and JH are based on the time-slotted 
round-trip carrier synchronization approach. 

After synchronization, all collaborative nodes transmit the signal coherently 

tr = Zke^^', r G M\ (7) 
Then the received signal at angle (p can be given as 

g{cf>) = Y.'^Y. y^rarke^'^e-^'^^^^ + w (8) 

where w ~ CA/'(0, cr^) is the additive white Gaussian noise (AWGN) at the direction 0. Note 
that the white noise is the same at all angles and, therefore, it disturbs the CB beampattem 
^ equally in all directions. The received noise power ct^ at BSs/APs can be measured in the 
absence of data transmission and, therefore, is assumed to be known at each BS/AP. 
The received signal at the intended BS/AP dk* can be written as 

9k* = gi^k*) = Zk* ^ A/^flrfc- + ^Zk^ A/^a^fc-e"^(^' + w 
where xi'^ '^^ = 71 ^^e~^^^^ "^^■'j, yi'' '^^ = xje"-'^^'^ -^r)| and TZ{-} and X{-} represent the 

(k* k) (k* k) 

real and the imaginary parts of a complex number, respectively. Note that Xr ' and yr ' are 
random variables. It can be further shown (see Appendix A) that u G |xr^ yr^ ''^■'j has mean 
m^(k*.k) = m^(k*,k) =mu = E{u} = Q and variance 0"^^.,^) = cr'^^^k-'.k) =(^l = E {u^} = 0.5. The 
first term in ^ is the signal received at the the BS/AP dk* from the desired set of collaborative 
nodes A^'^*, while the second term represents the interference caused by other sets of nodes 
J^'^, k ^ k* where Ai'^ fl Ai^ = 0, k ^ n. Assumed that each node in the network utilizes 
the same amount of power for each CB transmission, i.e., Pr = P, the received signal ^ at the 
BS/AP dk" can be rewritten as 
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ark* ^\fP^Zk ^ ark* {xf'^'^ -jyf'^^) + w. 



(10) 



k^k* r£M^ 



III. SiDELOBE Control via Node Selection 



In dense WSNs, each source node is surrounded by many candidate collaborative nodes in 
its coverage area. Although it has been shown earlier that the sidelobe levels of the average 
beampattern decrease inverse proportionally and uniformly over all directions with the increase 
of the number of collaborative nodes [|3l, flU, the sidelobe levels at particular directions of interest 
(unintended BSs/APs) in the sample beampattern are totaly random and can be unacceptably high 
if the number of collaborative nodes is not very large. At the same time, the randomness of the 
node locations provides additional degrees of freedom for controlling the beampattern sidelobes. 
Indeed, the sidelobes corresponding to different sets of collaborative nodes are different flSl, BU. 

The problem of high sidelobe levels at specific directions for the average beampattern has been 
briefly discussed in [|26l . It is suggested there to use only the sensor nodes placed in multiple 
concentric rings instead of using all nodes in the disk of the coverage area. However, a narrower 
ring with larger radiuses results in the average beampattern with smaller mainlobe width and 
leads to larger sidelobe peak levels at other uncontrolled directions than the conventional CB. 
Moreover, only the average beampattern behavior is considered in [|26il . while it is the sample 
beampattern behavior that is of real importance for the sidelobe control in WSNs. 

Exploiting the randomness of the node locations in WSNs, we introduce and study in this 
section a node selection algorithm for sidelobe control of the CB sample beampattern as a method 
for interference reduction. A sidelobe control approach based on node selection is suitable for 
WSN applications because it allows to avoid complex central beamforming weight design and 
corresponding additional communications. 

It is required to select a subset of collaborative nodes from the candidate nodes in the coverage 
area of the source node. Let M'' be a set of collaborative nodes to be selected from Ai'', i.e., 
M'^ C Ai'^, to beamform data symbols to dk- Note that in accordance with [3] and [4], the 
mainlobe of the beampattern is stable and does not change for different subsets of Ai'' as long 
as the the size of the coverage area does not change and the WSN is sufficiently dense, i.e., each 
cluster of the WSN consists of a sufficiently large number of sensor nodes. Also note that the 
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set M'' can be updated any time when the channel conditions or network configuration change. 
Alternatively, it can be updated periodically to balance power consumptions among nodes. The 
meaningful objective for node selection is to achieve a beampattern with low level sidelobes 
toward the unintended BS/AP directions. 

Toward this end, we develop a low-complexity distributed node selection algorithm, which 
guarantees that the sidelobe levels toward unintended direction(s) are below a certain prescribed 
value(s) as long as the WSN is sufficiently dense. An algorithm utilizes only the knowledge of 
the received interference power to noise ratio (INR), denoted as 77, at the unintended destinations 
and requires only low-rate (essentiality, one-bit) feedback from the unintended BSs/APs at each 
trial. Although our node selection strategy does not guarantee the optimum result of centralized 
beamforming strategies (which require global CSI and, therefore, a very significant data overhead 
in the network), it has the following practically important advantages for applying in the WSNs 
context. 

(i) It is very simple computationally and can be run in cheap sensor nodes without adding 
any computations. 

(ii) It has a distributed nature and, therefore, uses minimum control feedback from the unin- 
tended BSs/APs. 

In the following, we first describe the communication protocol and then give the details of 
the node selection algorithm. 

A. Communication Protocol 

Data transmission is organized in the following steps. 
Step 1: At the beginning of each time-slot, the source node listens to the control channels from 

BSs/APs and checks for an available BS/AP. 
Step 2: The source node Sk broadcasts its ID and the ID of the available BS/AP dk to the nodes in 
its coverage area Ai''. Note that all transmissions from the source node to the intended 
destination rf^ will be achieved through the CB using the nodes in the set of collaborative 
nodes M'' C M''. 

Step 3: The source node attempts to transmit its ID to the target BS/AP. If a collision occurs, 
i.e., if the target BS/AP receives transmissions form more than one collaborative sets of 
sensor nodes, the BS/AP approves only one set of nodes for data transmission. 
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Step 4: At the beginning of each following time-slot, the BS/AP broadcasts through the control 
channels the selected source ID in addition to one bit of information which indicates that 
the BS/AP is busy. In this case, other source nodes are not allowed to transmit data to the 
same BS/AP 

Step 5: Finally, only the predetermined subset of collaborative nodes assigned to the pair s^-dk 
continues to receive data, while other nodes go back to idle mode, that is, the actual data 
transmission takes a place for the pair Sk-dk- 

B. Node Selection Algorithm 

A set of collaborative nodes M'^ C A^'^ is assigned to each source-distention pair Sk-d^. 
To select such a collaborative set, the nodes can be tested one by one or a group of nodes by 
a group of nodes. The latter is, however, preferable since it can significantly reduce the data 
overhead in the system. Indeed, while testing one node or a group of nodes, we need to check 
if the corresponding CB beampattem sidelobe level reduces in the unintended direction(s) and 
then send the 'approve/reject' bit per one node in the fist case, or per a group of nodes in the 
second case. Therefore, if every group of nodes consists of a larger number of sensor nodes, 
less 'approve/reject' bits has to be sent in the system in total. 

Consider the source node Sk*, let the number of nodes in its coverage area is M, the number 
of collaborative nodes needed to be selected is < M, and the size of one group of nodes to 
be tested in each trial is L < A^. Using the selection principle highlighted above, the selection 
process can be organized in the following two steps. 
Step 1: Selection. Source node s^* initiates the node selection by broadcasting the select message 
to the nodes in its coverage area, namely the set Ai''* , and randomly selects a subset 
of L candidate nodes from Ai^* . 

The nodes can be assigned to the set C''* by using any of the following two methods. The 
first one is a centralized method in which the source node s^* is totally responsible for the 
node assignment. According to this method, every source node maintains a table of IDs of 
all candidate nodes in its coverage area and randomly assigns nodes to the set C''* . The 
source node s^* then broadcasts the IDs of the nodes assigned to the set C''* to inform 
them that they are selected for the test. The disadvantage of this method is that it is suitable 
only for small WSNs where each source node can keep records of all other nodes in its 
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coverage area. As a result, this method typically requires large data exchange between the 
source and the candidate nodes. 

Alternatively, in the second method, node assignment task is distributed among the source 
and collaborative nodes. In particular, if collaborative nodes receive the select message, each 
node starts a random delay using an internal timer. After the random delay, the candidate 
node responds by the ojfer message which contains the ID of this node. Then the source node 
responds by the approval message which requires only 1 bit of feedback. If a collision occurs 
and two collaborative nodes transmit the ojfer message at the same time, the source node 
responds by the approval message with a different bit value and the timers in both nodes 
start over a new random delay. The process repeats and the source node s^* keeps sending 
the select message until L candidate nodes are assigned and the set C^* is constructed. 
Step 2: Test. The set C^* transmits the test message that contains the intended BS/AP ID to the 
intended destination dk* using CB. While the intended destination dk* receives a predeter- 
mined signal power leveR the interference power levels at the unintended destination(s) dk, 
k ^ k* are random because of the random sidelobes of the CB beampattem. At this stage, 
all unintended BSs/APs with different IDs measure the received INR rj. If t] is higher than 
a predetermined threshold value ?7thr, the reject message is sent back to the candidate set 
CKln this case, the nodes in the candidate set C^* are all returned to the set of nodes J^^* 
and can be used in future trials. If no reject message is received from any of the unintended 
BSs/APs after a predetermined time, then the candidate set C''* is approved and each node 
from the candidate set C'^* stores the IDs of the source node Sk* and the destination dk*. 
Then the collaborative nodes assigned to serve the pair Sk*-dk* do not participate in future 
trials. In this way, we can avoid an overlap between sets of nodes serving different BSs/APs. 
In order to select collaborative nodes, the Selection and Test steps are repeated until N/L 
candidate sets Cf*, I = 1, 2, ... , N/L, are appro vedl^ Then the so obtained set of approved 
collaborative nodes is M'^* = IJ^ £f * . 

^The power level at the intended destination depends on the number of collaborative sensor nodes and the power of each of 
them. 

'it is assumed for simplicity that N/L is an integer number. If N/L is not integer, it is still easy to adjust the size of the 
candidate set in the last trial of the algorithm only. For example, the size of the last candidate set can be chosen to 
be equal to the reminder of N/L. 
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Node selection algorithm 
Initial values: 

A'' and L are predetermined at the Source Node Sfe*. 

rythr is predetermined at the unintended Destinations d^, fc = {0, 1, . . . , D}. 

1: At Sfe.: (Counter I <— 1). 
2: If (Counter / < f ), 



3: Then: { sw broadcasts the select message. 

4: A candidate set £f is constructed. 

5: Using CB, the nodes in £f transmit the test message.} 

6: Otherwise: {Go to 12.} 

7: At any dk, k ^ k": If ( The received INR r) > r^ihr), 

8: Then { sends the reject message to £f .} 

9: Else { No reject message is received. 

10: >Cf is approved and the corresponding nodes store the IDs of Sk* and dk*. 

11: At Sfe. : (Counter / •(- Counter I + 1). Go to 2.} 



12: Sfe. broadcasts the end message. 

TABLE I 

Table 1 : Node selection algorithm for CB sidelobe control. 



Once nodes are selected, i.e., M^* is constructed, the source node Sk* broadcasts the end 
message and no more candidate sets C^* is constructed. 
The pseudocode of the node selection algorithm is given in Table IIII-B[ 

IV. Performance Analysis 

In this section, we analyze the proposed node selection algorithm in terms of (i) the average 
number of trials required for selecting a set of collaborative nodes which guarantee low interfer- 
ence level for unintended BSs/APs and (ii) the complementary cumulative distribution function 
(CCDF) of the INR rj. The first characteristic allows to estimate the average run time of the 
algorithm, while the second characteristic is needed to estimate the achievable interference levels 
versus the corresponding interference threshold values. 
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A. The Average Number of Trials 

Note that the signal-to-noise ratio (SNR), denoted as 7, received through the link Sk*-dk* 
at the intended BS/AP should be above a certain level which guarantees the correct detection 
with high probability. This SNR at the intended BS/AP must be guaranteed regardless of the 
number of collaborative nodes participating in the CB. Therefore, we assume in our analysis 
that the total transmit power budget for each tested candidate set of nodes is kept the same in 
each trial of CB transmission. In particular, if the SNR at the intended BS/AP is required to 
be 101og]^o(7) then in the selection process, the power per one sensor node in £f* has to 
be set as P = 0-^7//^ < P™^"", where P™^^ is the maximum available power at the node. Also 
note that the power consumed for running the node selection algorithm is, in fact, proportional 
to the average number of trials in the algorithm. Therefore, it is preferable to construct a set of 
collaborative nodes with less number of trials. 

In order to derive the average number of trials for the node selection algorithm, we, first, 
need to find the probability that a candidate set of nodes Cf is approved as part of the set of 
collaborative nodes M^* . This probability is the same as the probability that the set £f * generates 
an acceptable interference at the unintended BSs/APs. Since we assumed that only one set of 
collaborative nodes M^* is constructed at a time, there is no interference present from other 
collaborative sets. Therefore, the interference power received at the unintended BS/AP dk from 
the tested candidate set of nodes £f* which targets the intended BS/AP df can be written as 

where Zk*, ark, x^r and yr^ ''^^ are defined in Section II. 
Equivalently, (fTTI) can be rewritten as 

r&Cf 

= * (/?E 4" " - .y^E y'i''A ('2) 

V re£f* reef J 
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where x'^^ ''^^ = arkxi'' ''^^ and y 



'^^ '''^ = arkui'' '^'^ with mean and variance given as 



mi 



1 1 

a-rk U 



2 2 



(13) 
(14) 



Let us introduce the notations 



{k*,k) ^ 



r{k* ,k) A_ 



X ^^g£fc* y'l^^ '^^ Then X\'^ and y^*-" can be approximated by Gaussian random variables 



^{k*,k) 



r{k*,k) 



||3l. [[TT| with mean 
finally rewritten as 



my = mi = and variance = cry = 7(T^(Ti. Thus, (fT2)) can be 



(fc*,fc) 



(15) 



Using (fTS]) and the fact that \zk*\'^ = 1, the received interference power at the unintended 
BS/AP from the candidate set of nodes £f * can be expressed as 



\i{^,\cf)f = (x^ 



-{k*,k) 



Y, 



(k*,k) 



(16) 



The probability that the candidate set of nodes £f * is approved to join the set of collaborative 
nodes M^* , i.e, the probability that the INR rj from £f* at the unintended BS/AP is lower 
than the threshold value r/thr, can be then found as 



Pr (?7 < ?7thr) 



Pr 



Vi^k\c^ 



Pr 



/ (-^{k*,k) 



< ?7thr 



1 — exp 



where the INR rj 



X. 



(k*,k) 



r{k*,k) 



+ yi 



(17) 



cr^ is exponentially distributed random vari- 



able with the probability density function (pdf) 



fiv 



2o\ 



1 2^ 



0, 



7] > 
7] <0 



(18) 



If D unintended BSs/APs are present in the neighborhood of the set of candidate collaborative 
nodes Cf, the probability that the INR from Cf* at any one of these unintended BSs/APs is lower 
than the threshold value ?7thr is given by (flTl) . Therefore, the probability that £f* is approved by 
all BSs/APs is the product of the probabilities that Cf is approved by each of the unintended 
BSs/APs, that is. 



P 



1 — exp 



D 



(19) 
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It can be seen from (fT9l) that p decreases if the threshold ?7thr decreases or the number D of 
neighboring BSs/APs increases. 

Using (fT9l ), a closed-form expression for the average number of trials required by the node 
selection algorithm can be derived. Note that the actual number of trials is a random variable 
which we will denote as T. In order to construct the set of collaborative nodes M^* , Tq = N/L 
trials must be successful among all T trials. Since the candidate nodes are selected randomly at 
each trial of the node selection algorithm, the algorithm itself can be viewed as a sequence of 
Bernoulli trials. Since To of these Bernoulli trials must be successful in order to construct A/''^,* 
the probability distribution of T in a sequence of Bernoulli trials is, in fact, negative binomial 
distribution, that is, 

Pr(T = t)= Q ~_\)/" (1 - pY"^" ■ (20) 

Using (|20l) . the average number of trials for the proposed node selection algorithm can be 
obtained as (see the details of the derivation in Appendix B) 

E{T} = ^ = ^. (21) 
p L ■ p 

It can be seen from (|2T]) that the average number of trials is proportional to the size of the 
set of collaborative nodes but it is inverse proportional to the size of the candidate set of 
nodes £f* and to the probability that the set £f* is approved to join the set N^* . Therefore, 
less number of trials is required in average for the proposed node selection algorithm if L is 
chosen to be large or N is small. Moreover, if the probability p, which, in turns, depends on the 
threshold value ?7thi- of the INR allowed at the unintended BSs/APs from the set £f* (see (fT9l)). 
is large, then less number of trials is required. 



B. The Complementary Cumulative Distribution Function of the Interference 

Let us assume that K{< D) collaborative sets are active and target different destination from 
dk", these sets are M^, k ^ k*, and their union is denoted hereafter as [JM'^^''* . It can be seen 
from (flOl) that the total interference collected at the destination dk* from these K collaborative 
sets is 

k^k* rGAf 
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where the power per one collaborative sensor node is P = cr'^'y/N because the SNR at the 
intended BS/AP must be 101ogio(7) dB. Using the fact that A/"^ = [J and multiplying and 
dividing the right hand side of (l22l) by the total interference at dk* can be expressed as 



k^k* 1=1 re£f 
N/L 



(^'{k-\k) _ ^ Ak*,k) 



(23) 



k^k* 1=1 



and f/'*''^) 

'{k*,k) 1 -^(fc*,fc) 



where and Y^'^ are zero mean truncated Gaussian distributed random variables corre- 



sponding to Xl ' and Yf- ' of (flSl) for only the approved candidate subsets. It can be shown 

is m 



that the marginal conditional probability density function of t//^ G <| X^"^ , Y^ 

2 



■(>(/£*, fc) 



f 



Uk',k) 



7] < 1], 



thr 



(k*,k) 



X, 



(fc*,fc) 



27rCT| 



X exp 



1 — exp 



+ 1^, 



,(fc*,fc) 



1-2 Q 



r{fc*,fc) 



\ 



2al 



|<^/ 



(24) 



where Q(x) = l/-\/27r J^°° exp(— m^/2) is the Q-function of the Gaussian distribution. 
Using (|23l) . the total INR at rffc. can be then expressed as 



at. 



L 



' N/L 



(fc*,/c) 



1=1 



' N/L 
«=1 



(fc*,fc) 




{k*,k) 




(k*,k) 



(25) 



Based on the central limit theorem, both the real and imaginary parts of the interference 

, ri-^N/Ly{k\k) 

'I ""^ \ N^ 2^1=1 ' are 



from each set of collaborative nodes, i.e., ^ / ^i=i 



yN/L ^{k-^k) 
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zero mean Gaussian distributed random variables with variance for both given as 




=1 




X 



(l-(l + /3)e-0 



(26) 



/(,|ir,a) = " -^Z: forK>0,,>0,a = ^. ill) 



where (3 = o-,^?7thr/2a|. 

Since the real and imaginary parts of the interference from each set of collaborative nodes are 
Gaussian distributed, the INR from each set of collaborative nodes is exponentially distributed for 
a given constant noise power c^. Therefore, the total INR t] collected at dk* from all collaborative 
sets is a sum of exponentially distributed random variables (see (|25l) ) and can be shown to be 
Erlang distributed, that is, 

a^{r])^~^ exp (—a?]) 

{K-iy. ' 2aj' 

Finally, using (l27l) . the CCDF of the INR, i.e., the INR distribution over the the sidelobes of 

the CB beampattern with node selection, can be expressed in closed-form as 

^Hv>Vo) = Y.^^ (28) 

where K{< D) is the number of active sensor node clusters in the neighborhood of the BS/AP 

dk*- 

V. Simulation results 

In order to demonstrate the advantages of the proposed node selection algorithm for the CB 
beampattern sidelobe control and verify the accuracy of the analytical expressions, we present 
also the following simulation results. 

A. Sample CB Beampattern 

In our first example, we study numerically the performance of the proposed node selection 
algorithm for the CB beampattern sidelobe control. Unless otherwise is specified, the nodes are 
assumed to be uniformly distributed over a disk with radius R = 2X. The total number of sensor 
nodes in the coverage area of the transmitting source node is M = 512 and the desired number 
of collaborative nodes to be selected is = 256. The power budget per all N collaborative 
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nodes equals to 20 dB. The size of a group of candidate sensor nodes L is taken to be equal to 
32 and the INR threshold value at the unintended BSs/APs is set to r/thr = 10 dB. 
We report the results for the following four different cases. 
Case 1: The intended BS/AP is located at the direction (^q = 0°. There are D = A unintended 
neighboring BSs/APs are present at the directions (fi = —160°, (p2 = —50°, (fs = 60°, and 
V?4 = 170°. 

In Fig. |2l we plot the sample beampattern corresponding to the CB with node selection 
and compare it to the sample beampattern corresponding to the CB without node selection 
and the average beampattern. The directions to the unintended BSs/APs are marked by 
symbol 'x'. It can be seen from the figure that the CB with node selection archives the 
lowest sidelobes in the directions of unintended BSs/APs, while the sidelobes of the CB 
without node selection are uncontrolled and high in the directions of unintended BSs/APs. 
Moreover, it can be seen that due to the CB coherent processing gain the peak of the 
mainlobe corresponds to 3 1 dB for all the beampatterns, while the power budget per all N 
collaborative nodes was set to 20 dB. The latter can be also predicted by the theoretically 
computed gain increase of log^^Q (see flU). 
Case 2: In this case, it is required that the data from different sensor nodes have to be sent to all 
4 BSs/APs simultaneously. Fig. [3] shows the CB beampatterns for the corresponding 4 CB 
clusters of 256 collaborative nodes selected from 512 nodes available in each cluster. Note 
that the sets of sensor nodes in all 4 CB clusters are different from each other and do not 
overlap. 

It can be seen from the figure that each beampattern has minimum interference at the 
direction of the mainlobes of the other beampatterns. Moreover, the mainlobes of the 
corresponding beampatterns all have the required mainlobes with a peak value of about 
31 dB. 

Case 3: In our third case, the neighboring BSs/APs are assumed to be located in the range (f) E 
[25° 45°] which is closed to the mainlobe direction of the intended BS/AP. The INR threshold 
value is set to r/thr = 10 dB. 

The beampattern of the CB with node selection and the average beampattern are shown in 
Fig. m It can be seen from the figure that the CB with node selection is able to achieve a 
beampattern with sufficiently low sidelobes over the whole range cf) G [25° 45°]. Note that 
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this case corresponds to the situation when the unintended AP is actually another cluster of 
sensor nodes distributed over space, which, therefore, cannot be viewed as a point in space. 
Case 4: In the last case, we assume that D = A neighboring untended Bs/APs are located at the 
angles corresponding to the peaks of the average beampattem, while the intended BS/AP 
is located at ^p^ = 0°. 

Fig. [5] shows the average beampattern and the beampattern of the CB with node selection. 
Note that the peaks of the average beampattern are located close to the mainlobe of the 
average beampattem. Therefore, the locations of the unintended BSs/APs are actually the 
worst locations in terms of the corresponding average interference levels. As it can be seen 
from the figure, using the node selection, we can achieve minimum interference levels at 
the directions of unintended BSs/APs even in this case. 
Summarizing, it can be concluded based on all these cases that there is an significant improvement 
achieved by using the node selection algorithm in reducing the sidelobe levels at the directions 
of unintended BSs/APs. 

B. Effect of The Algorithm Parameters 

The two parameters in the node selection algorithm are the INR threshold ?7thr and the size L 
of the candidate set of nodes CK 

In this example, it is assumed that the intended BS/AP is located at ipQ = 0° and there is one 
unintended neighboring BS/AP at the direction cpi = 65°. The noise power equals to cr^ = 0.05. 
The coverage area of the source node has a radius R = 5A. The INR threshold value changes 
in the range ?7thr = [—15 10] dB. The parameters of the Gaussin distribution corresponding to 
the lognormal distribution of the channel coefficients are m = and = 0.2. Monte Carlo 
simulations are carried over using 1000 runs to obtain average results. 

Fig. |6] demonstrates the effect of the threshold ?7thr on the average number of trials required 
to select the set of collaborative nodes M'' using the sets of candidate nodes of different sizes 
L G {16, 32, 64, 128}. It can be seen from the figure that the curves obtained using the closed- 
form expression (|2T)) for the number of trials are in good agreement with the simulation results. It 
can also be seen that by decreasing the threshold ?7thr, the number of trials increases. Moreover, 
the number of trials can be controlled using L. Indeed, as L increases, the number of trials 
decreases. It is important to note that because of the normalization factor in (fTT)) . the consumed 
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power at each trial is the same for different values of L and the total consumed power in the 
selection process is proportional to the number of trials. 

Fig. |7] shows the average interference level of the CB beampattem with node selection versus 
threshold ?/thi for different values of L G {16, 32, 64, 128}. It can be seen from the figure that the 
average interference level is proportional to the threshold ?7thr. Comparing Figs. [6] and |7] to each 
other, we can observe a tradeoff between the average number of trials required for node selection 
and the achieved average interference (sidelobe) level. It can be seen that lower interference level 
can be achieved by using smaller values of L at the expense of larger number of trials. 

In addition, it is worth noting that there is no limitations in the node selection algorithm on 
selecting the value of L as long as L < A^. However, the threshold ?7thr affects the average 
interference level and, therefore, depends on the sensitivity of the front-end receiver of the 
BS/AP. 

C. Ejfect of Number of Neighboring BSs/APs 

In this example, we study the effect of number of neighboring BSs/APs D to the performance 
of the node selection algorithm. Note that for simplicity we have always assumed that ?7thr is the 
same for all BSs/APs. Such set up does not restrict the generality of the proposed node selection 
algorithm since the selection is performed based on the accept/regect bit from the corresponding 
BS/AP, while ?7thr is used only in the BS/AP. Therefore, it is straightforward to use different 
threshold values at different BSs/APs, and no changes to the node selection algorithm have to 
be done. 

In Fig. m the average number of trials is plotted versus the threshold ?7thr for different values 
of £) G {1,2,3,4}. It can be seen from this figure that as D increases, the average number 
of trials of the node selection algorithm increases exponentially if the required INR threshold 
?7thr is low. Finally, it can be observed that the analytical and simulation results are in a good 
agreement with each other. 

D. The CCDF of the Beampattern Level 

In our last example, we investigate the CCDF of the beampattem level. 
Fig. [9] depicts the probability that the interference exceeds certain level, i.e., it shows the 
CCDF of interference for different values of the INR threshold r/thr G {~5, 0, 5, 10}. In addition. 
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Fig. [To] illustrates the CCDF of the interference for different numbers of active collaborative sets 
K E {1,2,3}. It can be seen from Fig. |9] that the CCDF of the interference increases as r/thr 
decreases. Moreover, as can be observed from Fig. \Wi the CCDF of the interference increases if 
K increases. The latter fact agrees with the intuition that for larger number of collaborative sets 
transmitting simultaneously, the overall received interference by all BSs/APs must be higher. 
The simulation results in both figures perfectly agree with our analytical results as well. 



Node selection is introduced for the CB sidelobe control in the context WSNs. A low-overhead 
and efficient node selection algorithm is developed and analyzed. In particular, the expressions 
for the average number of trials required for the proposed node selection algorithm and the 
CCDF of the beampattern sidelobe level of the CB with node selection are derived. The effect 
of the number of nodes selected at one trial to the algorithm performance is also investigated. It 
is shown that increasing the number of nodes selected at one trial reduces the number of trials 
at the expense of higher sidelobe levels. It is also shown that the CCDF of the beampattern level 
depends on the interference threshold value. From both the analytical and simulation results, 
we have seen that CB with node selection has perfect interference suppression capabilities as 
compared to the CB without node selection for which the beampattern sidelobes are uncontrolled 
and can cause significant interference to unintended BSs/APs. 

Appendix A: Derivation of the Mean and Variance of xi'"*''''' and yi''*''''' 



First, we find the probability distribution of n G <^Xr ,yr | and then find the corre- 
sponding mean and variance. Assuming that the angles 9'^ and 9'^* are uniform distributed in the 
interval [— 7r,7r], i.e., 9 ~ U[—n,n], it can be found that the difference A = 9^. — 9f has the 
following distribution 



VI. CONCLUSIONS 





(29) 



Using the equality u = 1Z {e^"^} 
in the interval [— 27r, 27r] as 



cos (A), we can find the roots of the equation u = cos (A) 



A 



1 



cos -^{u) 



A2 



TT + cos (u). 



(30) 
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Then the distribution of u can be found using the well known expression [|28l 

/(A) 



/(«) 



(31) 



Ai,A2 



lf(A)| 

where /'(A) is the first derivative of /(A). 

Substituting (|29l ) and (l30l) into (|3T1) . and using the facts that = |— sin(A)| and sin(cos"^(M)) 



a/I — n^, the distribution of n can be found as 

/(A) 



Ai,A2 



sin(A) 

27r — cos~^(m) ^ vr — (vr + cos~"^(m)) 



27r2 sin(cos~^(M)) 27r2 sin(7r + cos^^('u)) 
2ti 1 



27r2va^ 
The mean can be then easily found as 



TT 



rriu = E{u} 



uf{u)du 
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=rfn = —yl 
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(32) 



(33) 



And the variance can be found as 



f{u)du 



u 



1 ttvI 



zdu 



sin ^ix 



2n 



u 



2tt 



0.5. (34) 



Similar result can be also derived for m = I {^^^} = sin (A) in the same way. 



Appendix B: Derivation of (|2T]) 

Consider an infinite sequence of independent Bernoulli trials with probability of success p. Let 
Zi denotes the number of trials before the first successful trial. Then Zi is geometric distributed 
random variable Zi ~ Geom(p), that is, 

Vy{Zi = k) = {I - pf~^p, A; = 1,2,--- ,cx). 

The corresponding moment generating function (MGF) for (|35]) is 

oo 

MzM = E{e'^'} = ^e'\l-pf-^p 



(35) 



k=l 



1=0 ^ 



p)e^ 



(36) 
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Therefore, the average value of Zi can be found as 



(1 - (1 - p)e*)(pe*) - {{p - l)e*)(pe*) 



p)e 



t\2 



t=0 



(1- (l-p)e*)2 
p 1 

~2 ~ ~ 



t=0 



(37) 



p^ p 

Similarly, we can find the average number of trials between the first and second successful 
trials E{Z2}, second and third successful trials E{Z^}, and so on until Tq = N/ L successful 
trial. Since Zi, Z2, ■ ■ ■ , Zt^ are independent identically geometric distributed random variables, 
i.e., Zi ~ Geom(p), Ylfli^i negative geometric distributed, i.e., YlJli^i ~ NegBin(ro,p), 
with average 

To N 
p 



E\f:z. 

i=l 



ToE{Z,} 



L ■ p 



(38) 



where the fact that Zi, Z2, ■ ■ ■ , Z^^ are independent identically distributed is used again. 
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Idle Node 




Fig. 1. WSN model with multiple BSs/APs. 
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Fig. 2. Beampattem: The intended BS/AP is located at ipo — 0° and 4 neighboring BSs/APs at directions tpi — —160°, 
V?2 = -50°, (y32 = 60°, and ips = 170°. 
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Fig. 3. Beampattern: 4 simultaneous beampatterns with BSs/APs at directions tpi = —160°, (^2 = —50°, ip2 = 60°, and 
^3 = 170°. 
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Fig. 4. Beampattern: The neighboring BSs/APs are located in the range G [25° 45°] and rythr = 10 dB. 
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Fig. 5. Beampattern: The neighboring BSs/APs are at directions corresponding to the peaks of the average beampattern. 
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Fig. 6. Average number of trials E {T} versus threshold 77thr: M = 512, A'' — 256, fo — 0°, and ifii = 65°. 
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Fig. 9. The CCDF of the INR for different values of the threshold 7?,hr: M = 512, N = 256, L = 32, (po = 0°, and (pi = 65°. 
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Fig. 10. The CCDF of the INR for different values of K: M ^ 512, iV = 256, L = 32, ipo = 0°, and r?,hr = 10 dB. 
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